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Abstract. 

Near Feshbach resonances, n|ap ^ 1, systems of Bose and Fermi particles become 
strongly interacting/dense. In this unitary limit both bosons and fermions have very 
different properties than in a dilute gas, e.g., the energy per particle approach a 
value h^ii}!'^ jm times an universal many-body constant. Calculations based upon an 
approximate Jastrow wave function can quantitatively describe recent measurements 
of trapped Bose and Fermi atoms near Feshbach resonances. 

The pairing gap between attractive fermions also scales as A '--^ T?r?^^ jm near 
Feshbach resonances and is a large fraction of the Fermi energy - promising for 
observing BCS superfluidity in traps. Pairing undergoes several transitions depending 
on interaction strength and the number of particles in the trap and can also be 
compared to pairing in nuclei. 



1. Introduction 

Recent experiments probe systems of bosons PP and fermions 12113111131^] near Feshbach 
resonances by expansion and RF spectroscopy, where the interactions strengths and 
densities are large, n^/^|a| ^ 1. In the Hamiltonian for a bulk system of fermions or 
bosons 

^=E|^ + E-(r.-r,), (1) 

i=\ i<j 

interacting through two-body potentials v{r), one can then no longer apply the dilute 
limit pseudo-potential approximation: v{r) = 47r^^a5^(r)/m, where a is the scattering 
length, 

For example, the energy per particle in the dilute limit is 

— = 27rn^—, nap^<l, 2) 
IS m 

for bosons and half that for the interaction energy of fermions in two spin states. It 
would lead to infinitely large positive and negative energies at Feshbach resonances, 
a ±oo. In stead the strongly or dense limit also known as the unitarity limit is 
encountered where energies (and pairing gaps) are predicted [Z1IH| to approach 

-— = constant x n n\a\ ^ 1 , (3) 

N m 
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where the universal constant is a fundamental many-body parameter that only depends 
on the spin and the number of spin states of the particle. Several calculations 
and measurements of this constant now exist as will discussed in detail below. On 
dimensional grounds the energy per particle is expected to scale with density as n^/^ 
independent of the scattering length \a\ — cxo. Thus both the Bose and Fermi interaction 
energies "fermionize" in the unitarity limit. 

It is implicitly assumed that both the scattering length and the interparticle spacing 
tq = (3/47m)^/^ are much larger than the range R of interaction which is usually the 
case for the cold atomic clouds. It is also the case for low density neutron matter 
since two neutrons are just unbound with ^5*0 scattering length a = — 18 fm much larger 
numerically than the typical i? ~ 1 fm range of nucleon-nucleon interactions. In contrast 
the neutron and proton are weakly bound as the deuterium atomic nucleus. 

One of the predictions for the novel scaling laws in the unitarity limit for fermions 
[Tj and bosons P were based on a calculation using the Jastrow ansatz for the two- 
body correlations also referred to as the lowest order constrained variation (LOCV) 
approximation [^j. This is a very useful model as it extends from both the dilute to the 
unitarity limit and analytically continues across Feshbach resonances a ±oo. 

Studies of the transition from a Fermi gas of attractive atoms to a molecular BEG 
conventionally contain both atom and molecule components interacting resonantly (see, 
e.g., [El HH 1121 Uni HH CHI)- The Jastrow approximation incorporates the strongly 
interacting limit including the transformation to molecules very differently. We shall 
argue that the transition from a Fermi gas to molecular BEG is smooth in the sense that 
two-body correlations gradually build up across the Feshbach resonance approaching the 
two-body wave function of molecules consisting of two bound Fermi atoms. The large 
size of the molecules near Feshbach resonances are therefore included in the two-body 
correlation function between atoms. The transition is continuous as a ^ ±oo and will 
be described within LOGV. 

We shall briefly outline the LOGV approximation and calculate the energy per 
particle in detail for both Bose and Fermi atoms in sections 2 and 3 respectively. We 
will show that the results compares well to recent data on Bosons [T] and Fermions 
[21 ini E] also near Feshbach resonances. In the remaining sections 4-1-5 we will discuss 
BGS pairing in atomic trap with attractive Fermions and draw a connection to pairing 
in nuclei and nuclear matter. 

2. Bosons 

The LOGV method was developed for strongly interacting correlated fluids as '^He, ^He 
and nuclear matter in [HI and has more recently been applied to kaon condensation ^Hl 
and strongly interacting fermions [7j and bosons As explained in these references 
the Jastrow wave function 

vl/,(ri,...,r;v) = n/(r.-r,), (4) 

i<j 
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incorporates essential two-body correlations and is a good approximation for cold 
dilute and dense bose systems. We shall extend this wave function to fermions in 
the next section. The pair correlation function /(r) can be determined variationally by 
minimizing the expectation value of the energy, E/N = {"^IHI"^) / (\E'|\E'), which may be 
calculated by Monte Carlo methods that are fairly well approximated by including only 
two-body clusters. The basic idea of this method (LOCV) is that for r < tq the Jastrow 
function /(r) approximately obeys the Schrodinger equation for a pair of particles 

— +t;(r) rf{r) = Xrf{r) , (5) 

m dr'^ J 

where the eigenvalue energy of two atoms A = 2E/N . To take into account many-body 
effects, which become important when r is comparable to tq, the conditions f{r>d) = l 
and f'{r = d) = are imposed at the healing distance d, which is determined self 
consistently from number conservation 

n / p{r)d^r = 1 . (6) 
Jo 

In the dilute limit /(r) ~ 1 and so d = tq whereas in the dense limit d = (2/3)^/'^ro. 
Generally ~ tq. 

When the range R of the two-body potential v{r) is small compared to both \a\ 
and vq, the boundary condition at short distances is given by the scattering length 
{rf)'/rf = —1/a at r = 0. Solving the Schrodinger equation gives a wave-function 
r/(r) oc sm[k{r — b)] for positive energies. The boundary conditions and normalization 
determine the phase kb, the energy and the healing length d. When positive, the energy 
E/N = h'^k'^/2m is calculated in terms of the scattering length (see [S| for details) 

a tan k — 1 

- = ^ (7) 

d 1 + K tan K ' 

where k = kd. 

There is also a negative energy solution as shown in Fig. 1. For negative scattering 
lengths such many-boson systems are unstable towards collapse whereas the many- 
fermion system is stable as will be discussed below. Furthermore, the negative energy 
solution for both bosons and fermions, when the scattering length cross over from large 
and negative to positive, corresponds to the molecular BEC as we shall see shortly. 
The solution for wave function to the Schrodinger Eq. (0) is r/(r) oc sinh[A;(r — b)] or 
r/(r) oc cosh [A; (r — b)]. The energy E/N = —fi'^k'^ /2m is determined by 

a tanh n — 1 , . 

d 1 — K tanh k 

The multiple solutions to Eqs. (jHI) and ((Zj) corresponding to 0,1,2,3... number of nodes 
in the (two-body) wave-function /(r) are shown in Fig. (Q). In the limit R —>■ nodes 
inside the two-body potential and deeply bound states are irrelevant. 

In the dilute limit \a\ <^ ro Eqs. (jZj) and (jHI) give the correct energy per particle, Eq. 
(0). By contrast, in the unitarity limit \a\ ^ tq, Eq. ^ reduces to Ktann = —1 with 
solutions Ki = 2.798386.., K2 = 6.1212.., etc., and asymptotically /t„ = nvr for integer n. 
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Figure 1. LOCV calculation of the energy per particle {2mrQ)E/N ~ zLk^ in a BEC 
vs. scattering length. Dashed line is the dilute limit result whereas the dotted line is 
half the molecule energy in vacuum E/N = ~hy2ma^. 



The negative energy solution to Eq. (jH} reduces in the unitarity hmit to Ktanh k, = 1 for 
a ^ To with solution kq = 1.1997... As the scattering length cross over from — oo to +oo 
the negative energy state is analytically continued towards the molecular bound state 
with E/N = —h'^/2ma^ as a ^ +0. Therefore the LOCV approximation to many-body 
systems correctly reproduces the energies in both the dilute limit and the molecular 
limit. In the unitarity limit it correctly predicts qualitatively the n^/^ scaling of particle 
energies of Eq. 0. Quantitative comparison to experimental data for cold bose and 
fermi atoms in the unitarity limits will be given below. 

The energy of a repulsive BEC gas as it approaches a Feshbach resonance, a +oo, 
is of special interest 

El /3\2/3 ^2^2 2/3 

= - \ = 13.33/i2 . (9) 

These energies can be compared to the JILA experiment P] where the scattering 
length is tuned near the Feshbach resonance by a magnetic field B as 

a(5) = a,, (l - — — ^ ) . (10) 

V -D — ^Feshbach/ 

Coherent atom- molecule oscillations are induced in the BEC at a frequency v{B). We 
take this frequency as the difference between the mean field energy of two atoms in the 
atomic state and the molecular state 
2E^-E_ h 



= — 



h N 2Txm(P 



(4 - ■ (11) 
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Figure 2. Frequency of the atomic-molecular transition in a ^^Rb BEC as function of 
magnetic field near the Feshbach resonance. The full curve is the in vacuum transition 
frequency for two atoms to a molecule ly = ti/2TTma^. The LOCV prediction of Eq. 

for the transition of two atoms to a molecule in the many-body system is shown 
with dashed curve. Data from JILA p. 



Here, the indices -|- and - refer to the (lowest) positive and negative energy branches 
of Fig. (1) respectively as function of a. The negative energy branch approaches the 
molecular energy, i.e. E^/N —h^/2ma'^ as a — > -|-0. The resulting frequency is shown 
in Fig. 0. 

At the Feshbach resonance the LOCV frequencies approaches a constant value which 
according to Eqs. (Illj) is u = fi{K\ + K'^)/2'nm(P = 2.51hn'^^^ /m for the transition of 
two atoms in the unitarity limit to a molecular state. For the ^^Rb BEC of average 
density n = 2 x 10~^^cm~^ this results in a frequency u = 5.9kHz at the Feshbach 
resonance as seen in Fig. Q. The calculated transition frequency is in nice agreement 
with the JILA data close to the resonance. The calculated frequencies have a minimum 
just above the resonance of order ~ 4kHz because the atomic energy decrease faster 
with detuning (i.e., decreasing scattering length) than the molecular energy. 

The absence of experimental data in the 155-156G region is due to overdamping 
of the atom-molecule oscillations. The imaginary part of the self-energies may in this 
region be so large near the Feshbach resonance that the quasi-particles are not well 
defined. 

The effective range has been neglected in our calculations because we have assumed 
that the range of interaction is small as compared to scattering lengths and interparticle 
distances. The effective range of the interaction is by definition equal to the range for a 
square well potential. However, when the two-body interaction has a long tail and the 



Bosons and Fermions near Feshbach resonances 



6 



scattering length is not large, then corrections to the scattering amplitude and particle 
energies from the effective range appear. This is the case in Fig. (j21) off the resonance 
where the scattering length becomes small. Around -B^157G the frequency is slightly 
increased improving the agreement with data. 

The depletion of the condensate can also be calculated in the LOCV approximation 
from /(r). The condensate fraction is jSj 



3 , 



— ( sm K — K cos K) — 1 



(12) 



l-nf [l-/(r)]'rf3 
n J \rQ ^ 

which predicts quenching of the BEG for The condensate fraction differs from 

the dilute limit fraction: no/n = 1 — (4/-\/3vr)(a/ro)'^''^, as well as from the hard-sphere 
potential and exact quantum Monte Carlo calculations. 

The three-boson problem cannot generally be described by the scattering length 
only but depends on a three-body parameter as well ^H] which is referred to as non- 
universality. It was shown in ^n]; however, that for three bosons in a trap the energy 
was universal in both the dilute and in the unitarity limit - but not between these 
limits. Thus we may expect that the many-boson problem is universal as well in these 
two limits. The three-fermion problem is simpler in the sense that the Pauli exclusion 
principle prohibit Efimov states and therefore no three-body parameter enters, and the 
three- and many-body problem is universal. 



3. Fermions 

In the unitarity limit, n\a\^ S> 1, a system of fermions is dense and/or strongly 
interacting. The energy per particle has been calculated within Galitskii resummation, 
the LOCV method and recently by fixed node Green's function Monte Carlo (FN- 
GFMC). They find that the interaction energy does not extrapolate to ±cxd as the 
dilute result of Eq. Q does, but approaches a universal constant times the kinetic or 
Fermi energy, as in Eq. (jS)) 

^ = Ek,n + Ei„t = lEp[l + P]. (13) 
JM 5 

Here, (3 = Eint/E^in is a universal many-body parameter in the unitarity limit, which 
only depends on the number of spin states s through the Pauli principle and the 
density n = skp/Gn. For two spin states (3 is displayed in Table 1 for the various 
calculations and recent experiments as will be described in the following subsections. 
Recent measurements |21IS1I11E1IE] oi f3 seem to confirm this unitarity limit near Feshbach 
resonances. 



3.1. Galitskii resummation 

The scattering amplitude, which in the dilute limit is simply / = —a, can in the 
interacting many-body system be expanded in powers of scattering length via Galitskii 
resummation of ladder diagrams [20] . It was shown in [7] that resummation of scattering 
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Table 1. Calculations and measurements of the ratio of interaction to kinetic energy 
(3 = Eint/Ekin in the dense/strongly interacting regime a — s- — oo for two spin states. 





balc./bxp. 


Keierence 


0.26±0.07 


Duke exp. 


PI 


0.3-0.4 


ENS exp. 




0.4-0.6 


JILA exp. 


i 


0.68±0.1 


Innsbruck exp. 


i 


0.67 


Galitskii approx. 


m 


0.54 


LOCV approx. 


13 


0.56±0.01 


FN-GFMC 


m 



amplitudes lead to a f ^ kp scaling of the real part of the scattering amplitude in the 
unitarity limit. Brueckner, Bethe and Goldstone j2l] pioneered such studies for nuclear 
matter and ^He, which are more complicated cases since the range of interactions, 
scattering lengths and repulsive cores all are comparable in magnitude. In our case 
the range R of interaction is small, kpR <^ 1, and therefore all particle-hole diagrams 
are negligible. This is opposite to electro-magnetic interactions where the strength of 
the interaction is small but the range of interactions is large and must be screened by 
Debye or Landau damping implicit in loop diagrams. In the Galitskii resummation all 
higher order particle-particle and hole-hole diagrams are found to contribute to the same 
order when kplal ^ 1. Due to the very restricted phase space such higher order terms 
are usually neglected as in standard Brueckner theory. Truncating the resummation to 
second order one finds for two spin states |^ that j3 = —175/27(11 — 21n2) = —0.67 in 
the unitarity limit. 

3.2. LOCV 

The unitarity limit can also be studied in the LOCV approximation. For fermions the 
many-body wave function must be anti-symmetric. To first approximation one includes 
a Slater wave function ($), which is the ground state of free fermions 

vl>^s(ri,...,r^)=$n/(r.-r,'), (14) 

i<j' 

We refer to Ref. [23] where also a BCS wave function including pairing is employed. 
Because $ insures that same spins are spatially anti-symmetric, the Jastrow wave 
function only applies to particles with different spins (indicated by the primes). 
Therefore, the number conservation also applies to different spin states only and Eq. © 
should be multiplied by a factor (s — l)/s on the left hand side which affects the healing 
length. For example, in the unitarity limit a — oo the healing length now approaches 
d = ro{2s/{s - 1)3)1/3 = (37r/(s - l)y/^kp\ 
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Otherwise the LOCV calculation for fermions follows the lines as for bosons and 
the interaction energy is Eint = /2md'^ as calculated above. In addition the kinetic 
energy (3/5)-E'ir appears due to the Slater ground state. The total energy becomes 

where the ± refers to the positive and negative energy states discussed above. The 
energy becomes 0| 

E 3„ h-'nl „ [3 ^s- 1^2/3 

(16) 



:Ep — ^ = Ej, 



N 5 2mrf2 



5 V 37r 

With kq = 1.1997.. we thus obtain in the unitarity limit (3 = 0, —0.54, —0.85, —1.12, 
for s = 1, 2, 3, 4, ... respectively. 

For the gas to be stable towards collapse the energy must be positive, i.e., 1 + /3 > 0. 
Therefore the LOCV approximation predicts that up to s < 3 spin states are stable in 
the unitarity limit whereas in the Galitskii approximation only s = 1, 2 are stable. The 
marginal case s = 3 may be studied with ^Li atoms since they have three spin states with 
broad and close lying Feshbach resonances [3]. The stability of two spin states towards 
collapse has recently been confirmed for a ^Li and ^"K gases near Feshbach resonances. 
It has long been known that neutron star matter [221 with two spin states likewise has 
positive energy at all densities whereas for symmetric nuclear matter with two spin and 
two isospin states, i.e. s = 4, the energy per particle is negative. Nuclear matter is 
therefore unstable towards collapse and subsequent implosion, spinodal decomposition 
and fragmentation at subnuclear densities |24j. Above nuclear saturation densities, 
kpR^l, short range repulsion stabilizes matter up to maximum masses of neutron stars 
~ 2.2Mq, where gravitation makes such stars unstable towards collapse. 



3.3. Fixed-node Green's function Monte Carlo (FN-GFMC) 

In Ref. [25J the energy of a finite number of Fermions in a box has been calculated 
within FN-GFMC at a Feshbach resonance. Since the gas is in a meta-stable state with 
typical ~ 40 lower lying molecular states, the number of nodes in the Jastrow two-body 
wave-function must be fixed in Monte Carlo calculations. At the Feshbach resonance 
a —oo they find that the energy of the system increase with the number of particles 
as E = (1 + (3)NEkin with (3 = -0.56 ± 0.01. This is an upper limit since FN-GFMC is 
also a variational calculation of the energy. 

Pairing favors an even number of particles and from the odd-even staggering of 
E{N) a pairing gap A ^ O.bAEp is extracted from the FN-GFMC calculations at the 
Feshbach resonance. 

When comparing the various calculations of (3 it should be noted that the FN- 
GFMC value f3 = —0.56 includes pairing energies whereas the LOCV and the Galitskii 
approximations did not. Excluding pairing the FN-GFMC leads to a higher energy 
state with j3 = —0.44 [SHI- In comparison the expected contribution {3/8) A"^ / Ep from 
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Figure 3. Ratio (3 of the interaction energy over the kinetic energy vs. magnetic field 
for ^Li. Predictions from Eq. (|15|l with k calculated within LOCV are show with full 
and dashed curves for negative and positive scattering length respectively. The data, 
density and scattering length vs. magnetic field near the Feshbach resonance at 855G 
is taken from 

pairing energies in the dilute limit would lead to an even larger correction if the pairing 
gap in the unitarity limit is inserted. 

3.4- Experiments 

Several experiments with trapped Fermi atoms have recently measured energies in the 
strongly interacting or dense limit near Feshbach resonances. The energy in the trap 
(excluding that from the harmonic oscillator potential) is E/N = {3/8)Epy/l + p where 
Ep = {"iNf/^huo is the Fermi energy in a trapped non-interacting gas. 

The Duke group [2j measured the energy of ^Li Fermi atoms near a Feshbach 
resonance from expansion energies. After correcting for thermal energies, they find 
(3 = —0.26 ± 0.01 close to the Feshbach resonance kpa ^ —7.4. 

Bourdell et al. at ENS [3] measure expansion energies for ^Li near a Feshbach 
resonance. They find (3 = —0.4 ±0.1 as a — >■ —oo as shown in Fig. Q. The situation 
is more complicated at the other side of the resonance for positive a, where the ratio 
increases up to /3 = 0.2 ±0.1 in agreement with LOCV predictions, but then drops to 
a plateau close to the resonance. The LOCV calculations compare well to ENS data as 
a —OO whereas /3(a ±oo) exceeds the data especially at the plateau. Whether this 
is due to molecule formation jSHj or finite temperature, which can have a strong effect 
on P j2Z|, should be investigated. 
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The MIT experiment [3] on ^Li and the JILA experiment on ^'^K j3] measure 
the transition frequency between three hyperspin states which confirm the unitarity 
hmitation to energies and frequencies. The JILA experiment is performed at two 
densities and does not find a plateau. In extracting the JILA value for f3 as given 
in Table 1, the reduction factor (s — 1) in Eq. (fTHjl is replaced by s because of the 
specific population of spin states in their RF experiment. 

With the discovery of a molecular BEC the Innsbruck group [H] has been able to 
measure the size of the atomic cloud, which scales with (1 + PY^'^, around the Feshbach 
resonance at very low temperatures. They find (3 = —0.68 ±0.1 compatible with the 
theoretical estimates. 



4. Pairing of Fermi atoms in harmonic oscillator traps 

BCS pairing and superfluidity is expected for trapped Fermi atoms with attractive 
interactions. As atomic gases have widely tunable number of particles, densities, 
interaction strengths, temperatures, spin states, and other parameters, they hold great 
promise for a more general understanding of pairing phenomena in solids, metallic 
clusters, grains, nuclei, neutron stars and quark matter. 

In a uniform dilute gas with sufficiently strong interactions or large number of 
particles the bulk limit is reached, and the pairing gap is at zero temperature 

TT 



(17) 



.2akF. 

In the unitarity limit, kp\a\^l, the pairing gap approaches a finite fraction of the Fermi 
energy [Zj just as the energy per particle. The FN-GFMC calculations discussed above 
P3] find that the odd-even staggering energy or pairing gap in bulk is A ~ O.bAEp which 
is quite close the value A ~ 0.49£'i;' obtained by extrapolating Eq. (fTTj) to a — — oo for 
two spin states. 

Near a Feshbach resonance the strongly interacting Fermi gas becomes unstable 
towards molecule formation and the critical temperature Tc = (7/7?) A ~ 0.567A 
for BCS superfluidity is expected to cross-over towards the slightly smaller critical 
temperature for forming a Bose-Einstein condensate of molecules j30| I31j. If we take 
A = OMEf from FN-GFMC both the BCS and BEC critical temperatures, ~ 0.5A, 
are therefore around or above the lowest temperatures reported for trapped Fermi atoms 
[3 El in El El- This bodes well for observing BCS superfluidity in atomic traps and 
establishing the unitarity scaling of pairing gaps. 

As cooling improves further and smaller pairing gaps may be found, new pairing 
regimes take over that are important in other systems as solids, nuclei, neutron stars, etc. 
We shall outline the various pairing regimes shown in Fig. Q vs. density interaction 
strength and relate it to BCS pairing in bulk. 

The various pairing regimes shown in Fig. can, except from the dense or 
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Figure 4. Diagram displaying the regimes for the various pairing mechanisms (see 
text) at zero temperature in h.o. traps vs. the number of particles N = n^p/i and the 
interaction strength a < 0. The dotted lines indicate the transitions between single- 
shell pairing A G, multi-level, single-level, and multi-shell pairing. At the dashed 
line determined by 2G\n{'-^np) — huj the pairing gap is A ~ huj, and it marks the 
transition from multi-shell pairing to bulk superfluidity Eq. 1)17(1 . The pairing gap is 
A = O.dAEp above the full line n|a|^ > 1, which separates the dilute from the dense 
gas. From P2I 



unitarity limit, be calculated from in the dilute limit Hamiltonian 

^ = E f + Irnu^'^l) + 4-^^£ E ^^(r. " r,) , (18) 

i=l \ / i<j 

for atoms in a spherical harmonic oscillator (h.o.) potential. When energies are 
measured in units of hu, and lengths in Qosc = \Jh/mu} only two parameters remain in 
the Hamiltonian, namely the number of particles and the interaction strength as plotted 
in Fig. (jU 

When the traps contain relatively few atoms that are weakly interacting (lower left 
corner of Fig. (jH)) the level spectra are highly degenerate due to the SU(3) symmetry of 
the spherical h.o. potential. Due to the Pauli principle the h.o. shells n = 0, 1, ni? are 
filled up. Only the top Fermi shell ni? ~ (3A^)^/^, where is the number of Fermi atoms, 
may be partially filled. The levels with angular momentum / = np, np — 2, ....,1 or are 
all degenerate in the weakly interacting limit due to SU(3) symmetry. Consequently, 
pairing is enhanced as the gap generally increases with the number of states, and leads 
to a supergap 



32^/2^ \a\ 

For stronger interactions multi-shell pairing also takes place whereas for more particles 
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the stronger mean field cause level splitting, which reduce pairing towards single level 
pairing \S3\ OH E2] and shows up as a distinct shell structure with h.o. magic numbers. 

Varying the number of particles and interaction strengths thus allows us to 
investigate a wide range of pairing mechanisms that apply to many other physical 
systems, e.g., nuclei and neutron stars as will be discussed in the following section. 



5. Pairing in nuclei and nuclear matter 

The nuclear mean field is often approximated by a simple h.o. form and the residual 
effective pairing interaction by a delta force in order to obtain some qualitative insight 
into single particle levels, pairing, collective motion, etc., (see, e.g., [SSI EE])- We can 
therefore compare pairing in nuclei to that in traps as investigated above, once the h.o. 
potential is adjusted to describe nuclei. 

As described in jS21 the anharmonic nuclear mean field is stronger and of opposite 
sign than that in atomic traps. It also contains a strong spin-orbit force which change 
the magic numbers from the h.o. shells. Correcting for that the quasi-particle energies 
and pairing gaps were calculated by solving the Bogoliubov-deGennes equations which 
can be expressed as gap equation. The pairing depends only on one parameter, namely 
the effective scattering length in the Hamiltonian a = — 0.41fm that was obtained from 
fitting the pairing gaps to the odd-even staggering binding energies of neutrons and 
protons in nuclei (see Fig. The shell structure and the average pairing is well 

reproduced. If nuclei were to be placed in the pairing phase diagram of Fig. (jH) 
they would lie in the transition region between the single-shell, multi-shell and single- 
level pairing regions. Because the multi-shell increase pairing and single-level decrease 
pairing the average gaps are to a good approximation given by the supergap of Eq. 

Adjusting the h.o. frequency and oscillator length to that in nuclei with constant 
central densities of n = 0.15fm~^ the supergap becomes 

A-G = ^^-^ (20) 
0.41fm AV3 • ^^^) 

The supergap predicts a A~^^^ scaling in good agreement with nuclear pairing data. 
It differs slightly from the standard A~^^'^ scaling in the Bethe-Weiszacker liquid-drop 
formula. 

The supergap is robust in the sense that it does not depend on the level spectra or 
other details, and thus allows us to extract the effective scattering length directly. 

For very large nuclei multi-shell pairing becomes increasingly important and pairing 
approaches that in bulk matter. With the effective nucleon coupling constant a = 
— 0.41fm extracted from pairing gaps in finite nuclei, we can estimate the ^Sq pairing of 
both neutrons and protons in nuclear matter from Eq. ()17|1 

A^l.lMeV. (21) 

Neutron star matter has a wide range of densities and is very asymmetric, Z/A ~ 
0.1. One can attempt to estimate of the pairing gaps as function of density from the gap 
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Figure 5. Neutron pairing A(^^(7V) vs. the number of neutrons. The experimental 
values have been averaged over isotopes |3Z|- The calculated gaps A and quasi- 
particle energies E are obtained from the gap equation (see text) with effective coupling 
strength a = —0.41 fm. The supergap G is shown with dashed line. From 

in bulk, Eq. ()17p. with a ~ —0.41 fm. However, the effective interaction a is density 
dependent. At higher densities we expect the effective interaction to become repulsive 
as is the case for the nuclear mean field at a few times nuclear saturation density. At 
lower densities the effective scattering length should approach that in vacuum which for 
neutron-neutron scattering is a{^So) ~ —18 fm. 

Neutron stars are so cold that neutron and proton superfluidity is expected in most 
of the surface layers. The depinning of superfiuid vortices can explain neutron star 
glitches in star quakes |42j . 

6. Summary and Outlook 

We list the main results discussed above and add a few topics for future investigation: 

• In the unitarity limit near Feshbach resonances, n\a\^ ^ 1, where bosons and 
fermions interact strongly or are dense, the particle energies (and pairing gaps) 
scale like h^n?^^ /m times an universal constant. 

• For ^^Rb bosons recent experiments |lj confirm the unitarity limit and measure a 
transition frequency v ~ 5kHz near the Feshbach resonance in agreement with the 
predictions from LOCV. Experiments at other densities are underway to check the 
"fermion" n^^^ dependence of the energy per particle. 

• The unitarity limit has also been confirmed for Fermi atoms near Feshbach 
resonances. Measurements of the universal parameter (3 = Eint/Ekin — —0.5 as 
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a — s> — oo are compatible with theoretical predictions. 

• The pairing gap in traps with attracting Fermi atoms is A ^ O.bEp near Feshbach 
resonances and the critical temperature ~ 0.5 A. Thus superfluidity may already 
have been achieved in recent experiments. 

• Cooling further will eventually make measurements of smaller pairing gaps possible. 
By varying the number of particles and interaction strengths the pairing in h.o. 
traps with attractive Fermi atoms undergoes several pairing phases. These exhibit 
a wide range of pairing mechanisms that apply to many other physical systems. 

• Pairing of neutrons and protons in nuclei is similar to super-pairing in atomic traps. 
The supergap scales as A ~ 5.5 MeV/A^^^ with mass number. Nuclei also extend 
into the multi-level pairing regime seen as a strong shell dependence with magic 
numbers. 

• Large nuclei approach the ^5*0 pairing bulk in nuclear matter A ~ l.lMeV for both 
neutrons and protons. Pairing will, however, vary with density and asymmetry in 
neutron star matter. Low neutron densities is also dominated by a large negative 
scattering length a(^SQ) — 18fm due to a Feshbach resonance in the NN channel. 

• Mixing fermionic with bosonic atoms allow sympathetic cooling EHl IHOl HHl 1^ 

and thus to study weaker pairing. Induced interactions between fermions and 
bosons generally enhance pairing [2^. 

• A systems of fermions with attractive interactions and two spin states will not 
collapse as bosanovae and supernovae, not even near a Feshbach resonance because 
the kinetic energy dominates, 1 + f3 > 0. Three spin states is marginal and may 
be studied in ^Li where three broad Feshbach resonances lie close in magnetic field. 
If three or more spin states can be tuned resonantly, their attractive interaction 
energy should dominate as a — ^ — oo and the systems collapse as "Ferminovae" . 

• Optical lattices in current experiments have few atoms in each local trap and we 
thus expect superpairing which favors the insulator vs. conductor state. 

Tabletop experiments at low temperatures at Feshbach resonances extend our 
studies of dilute systems to include also strongly interacting and dense Fermi and Bose 
systems with unitarity scaling. It will provide new insight into a new scaling region 
as well as BCS pairing in atomic traps. At certain interaction strengths and particle 
number the pairing mechanisms and superfluidity in h.o. traps with cold atoms are 
similar to that in nuclei and neutron star matter. 
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